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Abstract 



We study a noninteracting supersymmetric model in an expanding FRW spacetime. 
A soft supersymmetry breaking induces a nonzero contribution to the vacuum energy 
density. A short distance cutoff of the order of Planck length provides a scale for the 
vacuum energy density comparable with the observed cosmological constant. Assuming 
the presence of a dark energy substance in addition to the vacuum fluctuations of the 
field an effective equation of state is derived in a selfconsistent approach. The effective 
equation of state is sensitive to the choice of the cut-off but no fine tuning is needed. 

> 
On 

m : 1 Introduction 

It is generally accepted that the cosmological constant term which was introduced ad-hoc in 
the Einstein-Hilbert action is actually related to the vacuum energy density of matter fields. 
Observational evidence for an accelerating expansion [TJ[2j [3], implies that the vacuum energy 
density dominates the total energy density today. The vacuum energy density estimated in 
a simple quantum field theory is by about 120 orders of magnitude larger than the value 
required by astrophysical and cosmological observations [4] so that extreme fine tuning is 
needed in order to make a cancellation up to 120 decimal places. Theoretically, it is possible 
that the cosmological constant is precisely zero and the acceleration of the universe expansion 
is attributed to the so called dark energy (DE), a fluid with sufficiently negative pressure, such 
that its magnitude exceeds 1/3 of the energy density. Nevertheless, even if such a substance 
exists, it is extremely difficult to tune the vacuum energy to be exactly zero. Hence, the 
fine tuning problem perisits unless there exists a symmetry principle that forbids a nonzero 
vacuum energy. Such principle is indeed provided by supersymmetry [5]. In field theory with 
exact supersymmetry the contributions of fermions and bosons to vacuum energy precisely 
cancel [6]. However, the supersymmetry in real world is not exact. 

A nonzero cosmological constant implies the de Sitter symmetry group of spacetime rather 
than the Poincare group which is the spacetime symmetry group of an exact supersymme- 
try. Hence, the structure of de Sitter spacetime automatically breaks the supersymmetry. 
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Conversely, a low energy supersymmetry breaking could in principle generate a nonzero cos- 
mological constant of an acceptable magnitude. Unfortunately, the scale of supersymmetry 
breaking required by the particle physics phenomenology must be of the order of 1 TeV or 
larger implying a cosmological constant too large by about 60 orders of magnitude. Some 
non-supersymmetric models with equal number of boson and fermion degrees of freedom 
have been constructed [7] so that all the divergent contributions to the vacuum energy den- 
sity cancel and a small finite contribution can be made comparable with the observed value 
of the cosmological constant. 

In this paper we investigate the fate of vacuum energy when an unbroken supersymmetric 
model is embedded in spatially flat, homogeneous and isotropic spacetime. In addition, 
we assume the presence of a dark energy type of substance obeying the equation of state 
Pde = ^Pdej with w < 0. Unlike in flat spacetime, the vacuum energy density turns out to 
be nonzero depending on background metric. Hence, the expansion is caused by a combined 
effect of both DE and vacuum fluctuations of the supersymmetric field. Solving the Friedman 
equations selfconsistently we find the effective equation of state of DE. In particular, we find 
the conditions for which the effective expansion becomes of de Sitter type. The contribution 
of the supersymmetric field fluctuations is found to be of the same order of magnitude as 
DE and no fine tuning is needed. 

We do not claim that our model describes a realistic scenario but it is tempting to spec- 
ulate along the lines described in an earlier paper [8] where a naive model of supersymmetry 
in de Sitter spacetime has been considered. Our working assumption is that the universe 
today contains DE and no matter apart from fluctuations of a supersymmetric vacuum as 
a relict of symmetry breaking in the early universe. Since the global geometry is non flat, 
the lack of Poincare symmetry will lift the Fermi-Bose degeneracy and the energy density of 
vacuum fluctuations will be nonzero. This type of "soft" supersymmetry breaking is similar 
to the supersymmetry breaking at finite temperature where the Fermi-Bose degeneracy is 
lifted by quantum statistics ([9] and references therin). 

The remainder of the paper is organized as follows. In section [2] we introduce a super- 
symmetric model in an expanding FRW universe. The calculations and results are presented 
in section [31 In section H] we discuss the effective DE equation of state. Concluding remarks 
are given in section In appendix |X] we review the covariant regularization schemes of the 
vacuum expectation value of the energy momentum tensor in flat spacetime. 



Here we consider a noninteracting Wess-Zumino supersymmetric model with N species and 
calculate the energy density of vacuum fluctuations in de Sitter spacetime. In general, the 
supersymmetric Lagrangian L for N chiral superfields has the form [10] 
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where the index i distinguishes the various left chiral superfields and W(<&) denotes the 
superpotential for which we take 

W($) = ^^rmQ&i. (2) 

i 

Eliminating auxiliary fields by equations of motion the Lagrangian ([T]) may be recast in the 
form 

C = [WtWi - m-l&l 2 + -Va'd^t - -rmVi*?) , (3) 

i ^ ' 

where fa are the complex scalar and ^ the Majorana spinor fields. For simplicity, from now 
on we suppress the dependence on the species index i. 

Next we assume a curved background spacetime geometry with metric g^. Spinors in 
curved spacetime are conveniently treated using the so called vierbein formalism. The metric 
is decomposed as 

g,u(x) = Vab e\e b u ; g""(x) = tj^V, (4) 
where the set of coefficients e a M is called the vierbein and 

e/ = W"e\, (5) 

is the inverse of the vierbein. Obviously, 

(7 = det^ = -(dete%) 2 . (6) 

The action may be written as 

S = J d*x^g-(£ B + £ F ), (7) 

where Lb and Lp are the boson and fermion Lagrangians, respectively. The Lagrangian for 
a complex scalar field may be expressed as the sum of the Lagrangians for two real fields 

= 5 I> M >>:, - ™v 2 ) ■ (8) 



The fermion part is given by [12 
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i=i 



C F = l - ~ *;m7 m *) - (9) 

where 7 M are the curved spacetime gamma matrices 

r = e/7 a , (10) 

with ordinary Dirac gamma matrices denoted by r ) a . Variation of (J2J) with respect to \& 
yields the Dirac equation in curved spacetime 

iY^^-m^ = 0. (11) 
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The covariant derivatives of the spinor are defined as 



= *,„ - r^, (12) 



* ;M = * M + MV, (13) 

where 



with the spin connection [TTJ 



r, = ^/ 6 [7 a ,7 fe ], (14) 



w ; ft = ->iV(ev-r;/ A ). (is) 



In FRW metric the vierbein is diagonal and in spatially flat FRW spacetime takes a simple 
form 

e a M = diag(l, a, a, a). (16) 
where a = a(t) is the cosmological expansion scale. 

3 Calculation of the vacuum energy density and pres- 
sure 

A spatially flat FRW metric is given by 

ds 2 = dt 2 - aitfdx 2 . (17) 

It is convenient to work in the conformal frame with metric 

ds 2 = a(ri) 2 (dr] 2 - dx 2 ), (18) 

where the proper time t of the isotropic observers, or cosmic time, is related to the conformal 
time 77 as 

dt = a(rj)drj. (19) 

In order to calculate the energy density and pressure of the vacuum fluctuations we need 
the vacuum expectation value of the energy-momentum tensor. The energy-momentum 
tensor is derived from S as [12] 



r _ VabC fj, _ rpF , -rB ,. )n , 



B 



where the boson and fermion parts are derived from the respective scalar and spinor La- 
grangians 



T*=Y,d^ l d u ^- gtlu £B, (21) 



i=l 

i 



T lv = 1 (^T(mVV) - ^GmT^VO • (22) 



Owing to the assumed homogeneity and isotropy of spacetime the calculation of the density 
and pressure requires the Tq component and the trace T£. Specifically for the metric ([IB]) 
we obtain 

T B °o = «b = £ (^W) 2 + ^W) 2 + \m^ 2 ) , (23) 
^ = £ ( V) a + i W) 2 + 2mV t2 ) , (24) 

i=l ^ ' 



= h f = -— (^y^ - (a^)yv) + ^™/^, (25) 

tF " = i m ^- (26) 

Assuming a general perfect fluid form of the vacuum expectation value of T M „ 

< >= (p + p)u^u u - pg^, (27) 
the energy density and pressure of the vacuum fluctuations are given by 

p = u»u v < T^y >, (28) 



p =-(p-<T%>), (29) 



1 

where is the velocity of the fluid and < A > denotes the vacuum expectation value of an 
operator A. In particular, for vacuum energy we expect 



< ir >= PAS r, (30) 

in accord with Lorentz invariance. In this case we have 

Pa = -Pa- (31) 

With this equation of state we reproduce empty-space Einsteins equations with a cosmolog- 
ical constant equal to 

A = 87rGp A . (32) 

In the following sections we make the calculations in comoving coordinates. In comoving 
coordinates equations (128]) and (129]) simplify to 

Pvac =< T ° >, (33) 

Pvac = ^<T -T;>, (34) 
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3.1 Scalar fields 

Next we consider quantum scalar fields in a spatially flat FRW spacetime with metric (fl8|) . 
Each real scalar field operator is decomposed as 

p( V , x) = J2 a ^ (Xk(v)e rk "a k + Xk^Ye-^al) , (35) 
k 

in full analogy with the standard flat-spacetime expression (1HH1) considered in appendix lAl 
The function Xk satisfies the field equation 

Xk + {"m 2 a 2 + k 2 - a"/a)xk = 0, (36) 

where the prime ' denotes a derivative with respect to the conformal time rj. In massless 
case the exact solutions to this equation may easily be found [12]. In particular, in de Sitter 
spacetime a" /a = 1/r/ 2 , and one finds positive frequency solutions 

Xk = -jL=e-^ (l-f). (37) 



VWk V k V 

The operators a k associated with these solutions annihilate the adiabatic vacuum in the 
asymptotic past (Bunch-Davies vacuum) [T21 [TB] . 

If m solutions to flHBI) may be constructed by making use of the WKB ansatz 

XM = 7=4= e~^ naW ^ d \ (38) 

where the function W k is found by solving ( )36l) iteratively up to an arbitrary order in adiabatic 
expansion [TT]. For our purpose we need the solution up to the 2nd order only which reads 



W k = uj k + ^ 2 \ (39) 

where 

cu k = ^/m 2 + k 2 /a 2 . (40) 

The general expression for the second order term is [TT] 



(2) _ 3 1 a 2 3 1a 3 k 2 a 2 1 k 2 a 5 A; 4 a 2 
8ui k a 2 Au k a \a 2 uj\a 2 \a 2 uj\a 8o 4 w fe 5 a 2 ' 

where the overdot denotes a derivative with respect to the cosmic time t. Then, equation 
( 1391) may be written as 

^ -l t ($ + d i) [l+0(mV " l)] - (42) 

or, using ( TT9l) . as 

W k = u k - — ^ [1 + O(m 2 /io 2 k )] . (43) 
u k a° 
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We can calculate now the vacuum expectation value of the 0-0 component and the trace of 
the boson energy-momentum tensor. Using (|9"T|) and the commutation properties of a k and 
4, from (IH and (EI} with $35} we find 



< T B l >- 



V 



d 3 k 
(2^)3 



l\ '|2 , 2 2| I 



< T >- 



V 



d 3 k 



l\ l \2 2 2| |2 2 2| |2\ 

(IXfcl - a u k \Xk\ -am \x k \ ) ■ 



Using([33J) and (BSD with M3J we obtain 



d 3 fc 



a 3 y (27r) 3 a;/ s 



' 2 la' 2 la /2 m 2 1 , 



a' 2 a" 



a' of" 



-2 + C( 



-4\ 



id 



(44) 
(45) 

(46) 



The first term in square brackets is identical to the flat spacetime result. The second term 
is a quadratically divergent contribution due to a non flat geometry, the next two terms are 
logarithmically divergent, and the rest is finite. Similarly, with the help of (|34"|) we find the 
boson contribution to the pressure 



d 3 k 



a 3 J (2n) 3 ujk 



3a 2 
1 

+ - 









6 V a 4 




y 


a' 2 a" 


a'a"'\ 


i 


a 7 


a 6 ) 





1 



,'2 



m 



0{u k 



(47) 



3.2 Spinor fields 



Next we proceed to quantize the fermions. The Dirac equation in curved spacetime may be 
derived from fl9]). Specifically for a spatially flat FRW metric we obtain 



3 a 



2a 



«7° [d + -- ^ + i~j j djty - mty = 0. 



(48) 



Rescaling the Majorana fermion field ^ as 

= a - 3/2 i), (49) 
and introducing the conformal time we obtain for ip the usual flat spacetime Dirac equation 

Z7 d v ij} + iy'djij) — amip = 0, (50) 

with time dependent effective mass am. The quantization of ip is now straightforward [TBI [17]. 
The Majorana field ip ma y be decomposed as usual 



ifrfax) = (uksivV^hs + v ks {r])t 



l kx 



bl 



ks / ' 



(51) 



k,s 



where the spinor Uk s may be expressed as 

1 



Uks 



+_am(k) 
ak (k4>s 



(52) 
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Here, the two-spinors <j) s are the helicity eigenstates which may be chosen as 



(53) 



The spinor v ks is related to Uk s by charge conjugation 

v ks = ij°J 2 (uks) T - (54) 

The norm of the spinors may be easily calculated 

1 1 

UksUks = -VksVks = y(amQ - iQ')(am( k + iQ - — k 2 \( k \ 2 - (55) 

The mode functions (k satisfy the equation 

Ck + (m 2 a 2 + k 2 - ima')( k = 0. (56) 

In addition, the functions (k satisfy the condition [17] 

k 2 \( k \ 2 + (amQ - tC k ')(amCk + iQ = C x . (57) 

It may be easily verified that the left-hand side of this equation is a constant of motion of 
equation fl56l) . The constant C\ is fixed by the normalization of the spinors and by the initial 
conditions. A natural assumption is that at t — (rj — — 1/H, a — 1) the solution behaves as 
a plane wave (k = C 2 e~' tEkt , where E k = \/k 2 + m 2 . This gives Cfc(0) = C 2 , Cfc(0) = —iC 2 E k , 
and hence C\ = 2C\Ek{v(i + Ek). From (|5"5"|) and floT|) we obtain 

UksUks = -VksVks = y (Ci - 2/c 2 |C fc | 2 ), (58) 

which at t — reads 

m 

UksUks = -VksVks = ^itttt ■ (59) 

V Ek 

For C\ = 1 this coincides with the standard flat spacetime normalization |12j . 

In massless case the solutions to ( |56l) are plane waves. For m ^ two methods have 
been used to solve (!56|) for a general spatially flat FRW spacetime: a) expanding in negative 
powers of Ek and solving a recursive set of differential equations [16] b) using a WKB ansatz 
similar to ( 138|) and the adiabatic expansion [XT] . 

By making use of the decomposition ( 15"TT) and the standard anti-commuting properties of 
the creation and annihilation operators, the vacuum expectation value of the 0-0 component 
( 125]) and of the trace (T26]) of the fermion energy-momentum tensor may be written as 

< T F ° >= ^ ^Vksiam - ki)v ks , (60) 

< T F ^ >= ^ aVTLVksVks ■ (61) 

k.s 
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Evaluating the expression under the sum and replacing the sum with an integral as in (l9~Tj) 
we obtain 



1 



d 3 k 
(2^)3 



[ik 2 (CkC k ' - CkCk) ~ am] 



< T l " > 
- fj, " 



d 3 k 

W) 



(l-2A; 2 |a| 2 ). 



(62) 



(63) 



The expressions under the integral sign in f l62l) and ( 1631) were calculated by Baacke and 
Patzold [IE]. We quote their results for the divergent contributions: 



d 3 k 



^ ' 

" div ~ a 4 



— E k — 

d 3 k 



2E, 



l)m 2 ^ (a 2 



l) 2 m 4 



(2tt)' 



a 2 m 2 



Eh 



aa m 



8EI 



4 4 

am 



+ 



a' 2 m 2 



+ 



2 4 

am 



(64) 



(65) 



Note that the first three terms in square brackets in (l64"j) are identical to the first three terms 

-2 

J k 



in the expansion of auj k = ^JE^ + a 2 m 2 — m 2 in powers of E k 2 . Hence we can write 



P 



-< i*l >- 



d 3 k 



a 3 J (27r) 3 0Uk 



(66) 



The first term in square brackets is precisely the flat spacetime vacuum energy of the fermion 
field. The second term is a logarithmically divergent contribution due to the FRW geometry 
and the last term is finite and vanishes in the flat-spacetime limit a' — > 0. Note that, as 
opposed to bosons, there is no quadratic divergence of the type a' 2 /uk- 
Similarly, from ( 165]) we obtain 



< T Fil >-- 
- m 



d 3 k 



m 2 + ~^ + % 4 ) 

4 a 3 bj\ k 



(2n) 3 cu k 

and using (|34|) we find the fermion contribution to the pressure 



(67) 



V 



d 3 k 



(27r)3o; fc 



Ik 2 1 a' 2 m 2 1 a"m 2 



(68) 



3.3 Putting it all together 

Assembling the boson and fermion contributions, the final expressions for the vacuum energy 
density and pressure of each chiral supermultiplet are 



P B + P F 



d 3 k 



a 3 J (2ir) 3 u:k 
1 f 2 a' 2 a" 



la' 2 5a' 2 m 2 
2~^F + 8 a 1 " 



a'a'" 



a 
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2 + K 4 ) 



(69) 
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p 77" 

p = p + p 



d 3 k 



a 3 J (2ir) 3 u!k 



1 / a' 2 a"\ 1 / a' 2 a"\ m 2 
6 +24 )rf 



A- 



1 / a'V a'a'"\ 1 ^ 

+7 2 — r H -"2+ (0 

4 V a 7 a J u i 



(70) 



The dominant contributions in (169]) and (ITU]) come from the leading terms in square brackets 
which diverge quadratically. Note that these quadratically divergent terms are due to bosons; 
fermions only provide a cancellation of all divergent and finite terms in the respective flat 
spacetime contributions of bosons or fermions. 

To make the results finite we need to regularize the integrals. We will use a simple 3-dim 
momentum cutoff regularization (recently dubbed "brute force" cut-oof regularization |18j ) 
which, as shown in appendix rX] may be regarded as a covariant regularization in a preferred 
Lorentz frame defined by the DE fluid. 

The advantage of this approach is a clear physical meaning of the regularization scheme: 
one discards the part of the momentum integral over those momenta where a different, yet 
unknown physics should occur. In this scheme a preferred Lorentz frame is invoked which is 
natural in a cosmological context where a preferred reference frame exists: the frame fixed by 
the CMB background or large scale matter distribution. A similar standpoint was advocated 
by Maggiore [13] and Mangano [TJ]. Furthermore, as we have already demonstrated, a 
supersymmetry provides a cancellation of all flat spacetime contributions irrespective of the 
regularization method one uses. 

We change the integration variable to the physical momentum p = k/a and introduce a 
cutoff of the order of the Planck mass A cut ~ mp\. The leading terms yield 

N n' 2 f Acut NA 2 n' 2 

I p4K1 + 0(p- 2 ))=^^(1 + 0(A„>Ac U ,)), (71) 

p ^ is? ( 3 £ - 4) ^ + o(A » in A » i) ) • <72) 

where N is the number of chiral species. Clearly, we do not obtain the vacuum equation of 
state ( 131]) as may have been expected as a consequence of a regularization that assumes the 
existence of a preferred Lorentz frame. 

In order to estimate the cutoff we first neglect background DE and assume that the total 
energy density p is given by (ITT]) . If we compare the first Friedman equation with (ITT]) 
keeping the leading term on the righthand side we find that our cutoff should satisfy 

Acut = ^m P1 . (73) 

It is worthwhile to note that several approaches [131 HI [20] IZE] with substantially different 
underlying philosophy have led to results similar to ( ITT]) . In particular, Cohen, Kaplan, and 
Nelson [H] have employed a cosmological horizon radius Rh = 1/H as a long distance cutoff 
and derived an upper bound 

^v<|# (T4) 
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from a holographic principle. Here, Auv and L denote the ultraviolet and long distance 
cutoffs, respectively. Our result would saturate the holographic bound (!74j) if we identify 
a' /a 2 = l/L. 

The closest approach to ours is that of Maggiore [13] and Sloth [21] who present a similar 
calculation of zero-point energy using massless boson fields only. The main difference in 
[13] with respect to ours is that the cancellation of the quartic contributions was done by 
hand on the basis of the procedure used previously in the literature with the so-called ADM 
mass. In our model, the cancellation of all (not only quartically divergent) flat spacetime 
contributions is naturally provided by supersymmetry. Another difference is that our results 
( I7T1) and ( 1721) are sufficiently general to allow a self consistent approach. 

The above consideration gives only an estimate for the cutoff. In the next section we 
give a self consistent treatment of the supersymmetric vacuum fluctuations in the presence 
of DE. 



4 Effective equation of state 

Since there is no way to precisely determine the cutoff, it is convenient to introduce a free 
dimensionless cutoff parameter A of order A < 1 such that 




A cut = AW — mpi. (75) 



The factor 1/yN is introduced to make the result independent of the number of species. If 
we reinstate the cosmic time t, equations ( ITTj) and (172|) become 

P = A 8^' (76) 

, 1 fa 2 d\ 
P = A 8^U- 2 «J- (77) 

Obviously, the pressure is negative if a 2 < 2ad. E.g., for a de Sitter expansion we find 
a 2 = ad and p = — p/3. This case was considered by Maggiore [13] who concluded that 
the vacuum fluctuations cannot (at least in his approach) be interpreted as a part of the 
cosmological constant because in the second Friedman equation the accelerating effects of 
pressure are canceled by those from the density. We shall see shortly that this conclusion is 
slightly altered in a selfconsistent approach to the effective equation of state. 

In addition to vacuum fluctuations of matter fields, we assume existence of DE char- 
acterized by the equation of state pde = ^Pde- The Friedman equations then take the 
form 

—n — -it^Pde + A-r , (78) 
cr 3 cr 

a 471" . (a 2 d\ 

- = --G p DE + 3p DE - A — - - . 79 
a 3 \a z a ) 

Introducing the effective equation of state 

PcS = W eS p eS , (80) 
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(81) 



2 A 

w eS = w + -- — 



(82) 



equations (178"]) and ( 1791) may be recast in the standard FRW form 




(83) 



a 



— G{1 + 3w e ff)Peff • 



(84) 



a 



Three remarks are in order. First, it is clear from (18T|) why we have chosen the cutoff 
parameter A less than 1. Second, it follows from (|82|) that the contribution of the vacuum 
fluctuations to the effective equation of state is always positive and hence it goes against 
acceleration! The third remark concerns the Bianchi identity which would not be respected 
if the vacuum fluctuations were the only source of gravity in Einstein's equations. However, 
because of the additional contribution to the energy-momentum tensor coming from DE, 
it is not necessary to have both contributions separately conserved. Since the effective 
pressure and energy density satisfy Einstein's field equations (183]) and (1841 . the combined 
energy-momentum is conserved and therfore the Bianchi identity is respected. In this way 
an interaction between the vacuum fluctuations and DE is implicitly assumed in the spirit 
of the two component model of Grande, Sola and Stefancic. [22J. 

It is worthwhile to analyze interesting cosmological solutions to equations (183]) and (184]) 
depending on the nature of DE given by the equation of state pde = w/>de- 

1. Consider first the case when there is no DE, i.e., when pde = Pde = 0. In this case 
equations (1781) and (1791) admit only a trivial solution a =const. Clearly, if A = 1, 
equation (ITS"]) becomes a trivial identity and equation (179]) implies a = 0. If A =fi 1, 
equations (1781) and (1791) are satisfied if and only if a = 0. Therefore, a =const is the 
only solution to (ITS]) and (1791) for any choice of A. In other words, FRW spacetime 
cannot be generated by vacuum fluctuations alone in an empty background. 

2. Another interesting special case is DE represented by a cosmological constant, i.e., for 
the equation of state pde = — Pde- It follows from (I8"2~]) that an accelerated expansion 
(w c s < —1/3) is achieved for any value of the cutoff parameter in the range < A < 
1/2. This case has also been discussed in [T3] IT4"]. 

3. A more general case is obtained if we only require accelerating expansion, i.e., if the 
effective equation of state satisfies w e g < —1/3. Then equation ( 182]) implies that the 
range — 1 < w < —1/3 is compatible with < A < 1/2, whereas w < — 1 would imply 
A > 1/2. In the latter case the DE equation of state violates the dominant energy 
condition. The fluid of which the equation of state violates the dominant energy 
condition was dubbed phantom energy [231 El] and has recently become a popular 
alternative to quintessence and cosmological constant [25] . 
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4. In the last example, we require that the background be de Sitter , i.e., w e g = — 1. In 
other words the effective equation of state describes an effective cosmological constant. 
From (182]) we find 



31- A 

Hence, this case may be realized only for a fluid with w < —1, i.e., for the phantom 
energy. We see that in a selfconsistent approach, unlike in the example discussed in 
[13] , a de Sitter expansion can be achieved as a result of a combined effect of DE and 
vacuum fluctuations. 



5 Conclusion 

We have calculated the contribution of supersymmetric fields to vacuum energy in spatially 
flat, homogeneous and isotropic spacetime. In addition to supersymmetric fields we have 
assumed existence of a substance obeying the equation of state pde = wpr>E, with w < 0. 
Unlike in flat spacetime, the vacuum fluctuations turn out to be nonzero depending on 
background metric. Combining effects of both dark energy and vacuum fluctuations of the 
supersymmetric field in a selfconsistent way we have found the effective equation of state. 
In particular, we have found the conditions for which the effective expansion becomes of de 
Sitter type. The contribution of the supersymmetric field fluctuations is of the same order 
of magnitude as DE and no fine tuning is needed. 

We have found that if we impose a UV cutoff of the order mp\ the leading term in 
the energy density of vacuum fluctuations is of the order H 2 m^ h where H = a/a. In this 
way, if we identify the expansion parameter H with the Hubble parameter today, the model 
provides a phenomenologically acceptable value of the vacuum energy density. We have also 
found that a consistency with the Friedman equations implies that a natural cutoff must be 
inversely proportional to yN. A similar natural cutoff has been recently proposed in order 
to resolve the so called species problem of black- hole entropy [26] . 

A Covariant regularization of T^ v in flat spacetime 

To illustrate problems related to the field theoretical calculation of vacuum energy we review 
the well known results for the scalar field in flat spacetime [271 1281 [29] . Consider a single 
noninteracting real scalar field described by the Lagrangian 

£ = \rT<P#<Pr ~ \™\\ (86) 
with the corresponding energy-momentum tensor 

Tpv = dn(pd v (p - rj^uC. (87) 
The field operator is decomposed as 

f ) = E ~WW k [f lEkt ^^ + , (88) 
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where 

E k = \/m 2 + k 2 . (89) 

and at and a\ are the annihilation and creation operators, respectively, associated with the 
plane wave solutions with the standard commutation properties. 

k,4] = %*'- (9°) 

From ()86l)-(l88l) with ( 1901) and replacing the sum over momenta by an integral in the usual 

way 

(91) 

k 

we find the vacuum expectation value of T, 



1 r d 3 k 

2 J (2n) 3 E k 



< T /->= o / T^T^M^ ( 92 ) 



where k^ = (E k , k). The righthand side of (|92|) may be expressed in a manifestly covariant 
way 



< T„ v >= j ^L ktl kJ(k% - m 2 )9(k ) . (93) 

The delta function under the integral restricts the integration to the hypersurface defined by 

- m 2 = ; k > 0, (94) 

with the invariant measure d 3 k/Ek on the hypersurface. Performing the integral over ko in 
( 1931) one recovers ( |92l) . However, if one assumes the vacuum expectation value of T^ u to be 
of the form 

as dictated by Lorentz invariance of the vacuum, one encounters inconsistency since different 
results for p vac are obtained depending on which component of T^ u one calculates. For 
example, using T 00 one finds 

1 f d 3 k 

Pvac =< r 00 >= - y —E k . (96) 

On the other hand, using the trace one finds 

a =-T»=—( ^ (97) 

which does not agree with fl96l) . One must conclude that the assumption fl95l) is not com- 
patible with ( 193]) . The reason for this inconsistency is that the integrals in expressions fl92l) 
and (I9"3"|) are divergent and make sense only if they are regularized. 

One way to covariantly regularize f[9"2"j) or ( ]9~3"]) is to cut the hypersurface (J54]l by a spacelike 
hyperplane defined by 

= - VK 2 + m 2 = 0, (98) 
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where K is an arbitrary constant of dimension of mass and is a general future directed 
timelike unit vector, which may be parameterized as 



Ufj. = (cosh a, sinh a sin 9 cos 0, sinh a sin 9 sin <fr , sinh a cos 9) . (99) 

Clearly, the vector u p is normal to the hypersurface / = const because df /dk^ — u^. In 
this way one effectively introduces a preferred Lorentz frame defined by the vector u p as if 
the vacuum fluctuations are embedded in a homogeneous fluid moving with the velocity u^. 
The special form of the constant in ( |98l) is chosen for convenience. 

The hyperplane cuts the hypersurface f l94l) at a 2-dimensional intersection defined by 
( 1981) together with ( 1911) . This gives a quadratic equation the solutions of which define a 
2-dimensional closed surface as a boundary of the integration domain £ defined by 



VK 2 + m 2 - u p k^ > (100) 
together with (19^1) . Hence, the regularized expression for < T Mi/ > is given by 

1 f d 3 k 

or in a manifestly covariant form 

< T ^ >= f Sh^h5(k p k p - m 2 )9(u p k p )9(VK 2 + m 2 - u p k p ) . (102) 
Using a general perfect fluid form fT2"T|) , p and p are given by the invariant expressions 

»=\Lwk w - (l03) 

In comoving frame (a = 0) the integration domain E becomes a ball of radius K and we 
obtain 

^ r ->=5L^' (105) 

p =<r«>=- / , ^ fc 2 . (106) 

Hence, the described covariant regularization is equivalent to a simple 3-dim momentum 
cut-oof procedure. The integration yields 

K 4 m 2 K 2 IK 2 . . 

+ -r^r-^l n -^ + -> (107) 



167T 2 167T 2 647T 2 m 

1 if 4 lm 2 if 2 IK 2 . , 

P= 316^-3T6^ + 64^ ln ^ + -' (1 ° 8) 



where the ellipses denote the finite terms. 



15 



This result reveals two problems. The first one concerns the fine tuning. Assuming that 
the ordinary field theory is valid up to the scale of quantum gravity, i.e. the Planck scale, 
the leading term in ( 11071) yields 



P 




10 73 GeV 4 , 



(109) 



compared with the observed value 



Per 



10" 47 GeV 4 . 



(110) 



This huge discrepancy may be easily rectified in flat spacetime simply by subtracting all di- 
vergent contributions and redefining the vacuum to have its energy exactly zero. However, as 
soon as we demand that vacuum energy or cosmological constant is nonzero the calculations 
should be repeated in curved spacetime (e.g. de Sitter spacetime) and a simple subtraction 
of vacuum energy by fiat cannot be done. 

If, in addition to the vacuum fluctuations of the field, one assumes that there exists an 
independent cosmological constant term A, as a result one would find an effective vacuum 
energy 



In order to reproduce the observed value one needs a cancellation of the two terms on the 
right hand side up to 120 decimal places! The problem is actually much more severe as there 
are many contributions to vacuum energy from different fields with different interactions and 
all these contribution must somehow cancel to give the observed vacuum energy density. 

The second problem is related to the equation of state. Obviously, equations (11071) and 
(11081) do not reproduce the expected vacuum energy equation of state (13T1) . as required by 
Lorentz invariance. Instead we find p = p/3 for the quartic term, p = — p/3 for the quadratic 
term, and only the logarithmic term satisfies (I3ip . This violation of Lorentz invariance is 
not surprising since the adopted covariant regularization procedure assumes existence of a 
preferred Lorentz frame. 

In principle, it is possible to regularize the energy momentum tensor by imposing (130]) and 
ignoring the mentioned inconsistency of the derived covariant expression (I9"2"|) . Then, using 
the manifestly covariant form f[9"3"j) of the energy momentum tensor one can calculate the 
components using covariant regularization schemes which do not invoke a preferred Lorentz 
frame. For example, the dimensional regularization with the MS prescription gives [28] 




and one would conclude that a covariant regularization removes the Lorentz violating quartic 
and quadratic divergences and retains only the logarithmically divergent term which agrees 
with the logarithmic term of the 3-dim cutoff procedure in (11071) and (11081) . However, in the 
Pauli Villars regularization one finds [27] 



Peff = P + Pa • 






(112) 



Ppv = -Ppv = 



64vr 2 



1 




-K 4 + 2m 2 K 2 - m 4 ( In ^- + -) 



(113) 



16 



so in this covariant procedure the quartic and quadratic divergences are present with coeffi- 
cients different from those of the 3-dim cutoff procedure. Again, the logarithmic term agrees 
with that of (11071) and ( 11081) . Both dimensional and Pauli Villars regularizations have an 
unpleasant property that the leading term contribution yields p < 0. This property is un- 
physical since p =< Tq > should be positive for the scalar field as follows from ( 1231) . Ossola 
and Sirlin have argued [27] that the quartic term in f II 1 3 j) may be removed by demanding 
strict scale invariance in the limit m — > or by invoking the Feynman regulator. 

Two other Lorentz invariant regularization schemes were considered by Andrianov et al 
|29j : ^-function regularization and the UV cutoff regularization of the large wave-number 
field modes. It was concluded that the former method is not adequate in treating the 
cosmological constant problem as it redirects the problem from UV to IR region. The latter 
method with a suitable choice of the large wave-number cutoff reproduces the Pauli Villars 
regularization result f 1 1 1 3 j) . With the choice advocated in [29J one can get rid off the quartic 
term but then the coefficient of the quadratic term changes. 

We see from the above analysis that a covariant regularization is ambiguous although in 
all mentioned covariant methods the logarithmic term comes with the same coefficient as in 
the 3-dim cutoff procedure. With the exception of the dimensional regularization where the 
power low divergences are absent by definition, the quadratic term is allways present with a 
coefficient that depends on the regularization method. 
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